RKKY interaction and Kondo screening cloud for strongly correlated electrons 
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The RKKY law and the Kondo screening cloud around a magnetic impurity are investigated for correlated electrons in ID 
(Luttinger liquid). We find slow algebraic distance dependences, with a crossover between both types of behavior. Monte Carlo 
simulations have been developed to study this crossover. In the strong coupling regime, the Knight shift is shown to increase 
with distance due to correlations. 
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Since its discovery, the indirect Ruderman-Kittel- 
Kasuya-Yosida (RKKY) exchange interaction between 
localized magnetic impurities embedded in a host metal 
has played an important role in the theory of magnetism. 
The magnetic moment of one impurity scatters conduc- 
tion electrons, which are then seen by some other im- 
purity. This second-order process results in the 2k p- 
oscillatory RKKY interaction between different magnetic 
moments where kp is the Fermi momentum. For 

a lattice of magnetic impurities, this interaction favors 
magnetically ordered phases determined by the lattice 
geometry. The RKKY interaction is a basic ingredient 
for many phenomena in strongly correlated systems, e.g., 
magnetic impurities in quantum wires |^], normal-state 
magnetism in high-temperature superconductors or 
magnetic ordering in heavy fermion materials 

For uncorrelated conduction electrons, the RKKY in- 
teraction in d dimensions is ^ cos[2kpx]/x'^ , where x is 
the distance between the localized moments . Within 
the random-phase approximation or Fermi liquid theory, 
this law is not expected to change qualitatively ||^. At 
the same time, however, it is known that Coulomb in- 
teractions can modify spin- or charge-density correlation 
exponents in ID p-[l4|. In this Letter, we show that the 
RKKY interaction indeed exhibits only a slow algebraic 
decay cos[2kpx]/x^'' , with an interaction-dependent 
exponent gc < 1- To examine how the RKKY law is af- 
fected by the Kondo effect in such a strongly correlated 
system, we study the magnetic screening cloud around 
a single Kondo impurity. We predict the asymptotic 
behavior far away from the impurity and qualitatively 
discuss how RKKY relates to Kondo screening physics 
in a Luttinger liquid. The crossover between these two 
regimes has been analyzed by Monte Carlo (MC) simu- 
lation. 

To describe the low-energy properties of correlated ID 
conduction electrons, we employ the bosonization tech- 
nique [ppiol] . The spin-i electron field operator can equiv- 
alently be expressed in terms of spin and charge boson 
fields, which obey the algebra (we put h = 1) 
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where i,j denote the charge (c) or spin (s) degrees of 
freedom. The canonical momentum for the 9i phase field 
is Ili{x) = dx<j>i{x). Written in terms of the boson fields, 
the right- or left-moving component (p — ±) of the elec- 
tron annihilation operator for spin a = ± is 
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where a = vp/ujc is a short-distance cutoff (wc is the 
bandwidth cutoff, say, the Fermi energy, and vp is the 
Fermi velocity). The unitary zero- mode operators rjpa 
annihilate a particle from branch pa and ensure that an- 
ticommutation relations hold between operators with dif- 
ferent pa (lO|. In contrast to models without spin flips, 
they have to be considered explicitly here to account for 
all minus signs. 

The archetypical low-energy theory for correlated elec- 
trons in ID is the Luttinger liquid model pO| ] , which uni- 
fies the low-temperature physics of microscopic lattice 
models for strongly correlated fermions. There are only 
two relevant interaction constants gc and gs- The charge 
interaction constant is gc ~ [1 + 2U/ttvp]~^^'^ < 1, where 
U is the forward-scattering amplitude of the screened 
Coulomb interaction potential. The Luttinger liquid 
model assumes that one is away from half-filling so that 
Umklapp scattering is not present. In addition, electron- 
electron backscattering processes are neglected, albeit 
one can incorporate them by a renormalization of the 
interaction constants or by a perturbative renormaliza- 
tion group (RG) scheme, where the fixed-point value is 
given by 5s = 1 Therefore, we will put = 1 in the 
following to respect the underlying SU(2) spin symmetry 
of the electrons. The Hamiltonian of the clean system is 
then given by 
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where Vj = vp/gj is the velocity of charge or spin density 
waves for the case of full Galilean translational invariance 
considered here jl^ . 

Let us now add a spin-^ magnetic impurity at a; = 0. 
We use the standard contact term with the conduction 
electrons, Hj = Js(0)5 where J is the direct 

exchange coupling, S the impurity spin operator, and 
s{x) the spin density operator, which from Eq. (|^) reads 
in bosonized form 

Szix) = ^-L H cos[2/ci?x + \/27r 6'c(a;)] cos[V27r ^^(a;)] 

VSvr TTfl 

s±{x) = — exp[±\/27r i0s(x)]|±icrj, cos[\/27r0s(a;)] + 

+ cos[2kFX + V2Tr0cix)]j , (4) 

with s± = Sx i iSy. Here we have used that the rjpa 
show up only as bilinear forms, for which a convenient 
representation can be found in terms of Pauli matrices. 



local screening properties induced in a Luttinger liquid 
by the presence of a single impurity. A related quantity 
of direct experimental relevance is the local susceptibil- 
ity x{^) — d{sz{x)) /9-B, which was recently reconsidered 
for the uncorrelated case and is proportional to the 
Knight shift. Linear response theory gives {j3 — 1/kBT) 
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This replacement gives the correct sign for all possible 
products of rjjj^rip,^, pairs allowing for a nonvanishing 
contribution. Therefore the chosen representation is suf- 
ficient for the calculation of correlation functions involv- 
ing only spin or charge densities [T^ . 

To eliminate an explicit dependence of the interaction 
part Hi on the (j)s field, we perform a standard unitary 
transformation jl^, U = exp[\/27ri(^s(0)S'2], such that 
our final Hamiltonian reads 



UHU-^ = Ho 
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SzOxOsiO) + — (axSx cos[V2^0c] 
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X{x) = pC{x) + p / dx'{sz{x)sz{x')) 
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where the second part does not contribute in the pertur- 
bative RKKY regime. 

Since the slowly varying part of C'{x) leads only to sub- 
leading terms for x ^ a we restrict ourselves to 
the 2k p part in the following. After the unitary trans- 
formation, we obtain 

C{x) = cos[2fcf a:] 3i^[y2^g^(^)] coa[V2^9,{x)]) . (8) 
zna 

As can be seen from Eq. (^), there is no phase shift in 
the cos[2kFx] term. This is clear since Eq. (|^) docs not 
include elastic potential scattering by the impurity. 

The standard treatment of the RKKY interaction Q| 
corresponds to a calculation of the correlation function 
C{x) by lowest-order perturbation theory in the exchange 
coupling J. The finite-temperature result is {x 3> a) 
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For X <C XT, where xt = vp/ksT denotes the thermal 
lengthscale, this yields the RKKY law 



+ UySy cos[V27r6ls] + a^Sz cos[V27r6's] cos[V27re'c] , (6) 

/ x— 

where J = J — 2ttvf. The four interaction terms include 
two forward and two backward scattering terms with or 
without spin flip, respectively. Backward scattering 
coa[^/2ndc]) is responsible for RKKY oscillations, while 
Kondo screening arises due to spin flip terms S^/y)- 

Our subsequent discussion is based on the correlation 
function C{x) = {sz{x)Sz), where the brackets indicate 
a thermal average. This function describes the spatial 
correlation of the electron spin density with the impu- 
rity spin. Another impurity spin located at x would see 
this correlation, and lowest-order perturbation theory in 
J constitutes an exact derivation of the RKKY law 
While a quantitative discussion of the complicated inter- 
play between the RKKY interaction and Kondo screening 
behaviors requires a study of higher-order terms in the 
corresponding two-impurity model j7|J|] , the main quali- 
tative features of this interplay can already be extracted 
from C{x) [0. In this work, we therefore focus on the 
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while for x ^ xt, an exponential decay on the scale xt 
is obtained. It is obvious that spin-fiip events do not 
contribute to the perturbative result (|9). Therefore, to 
lowest order in J, one could just as well consider a static 
impurity, or, equivalently, a point-like magnetic field act- 
ing at X = 0. The presence of such a field induces 2kF- 
periodic oscillations in the spin density of the electrons, 
which are then responsible for the RKKY interaction. 
Thus the range function describing the decay of the 
RKKY oscillation amplitude displays only a slow alge- 
braic ^ x^^" law in the low-tcmperature regime x <^ xt- 
In the noninteracting case, gc — 1, the usual x~^ decay is 
recovered. This modification of the range function might 
come as a surprise, since the Coulomb interaction does 
not couple to spin densities. The slower decay is a many- 
body effect induced by the presence of correlations. 

Starting from order on, spin fiips contribute and 
it becomes mandatory to treat the dynamics of the im- 
purity spin. The most important aspect of the impurity 
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dynamics is the Kondo effect, leading to a screening of 
the impurity spin by the Luttinger Uquid spin density be- 
low the Kondo temperature jl^, Tk ^ J'^K^-Sa) _ Kondo 
screening of the impurity becomes important for strong 
couplings J or at low temperatures. For instance, the 
second-order contribution to C(x) dX x <^ xt is 
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The logarithmic corrections over Eq. ( p^ are typical for 
the Kondo effect and indicate that we are dealing with a 
nonperturbative problem. 

To study the crossover from the RKKY law to the 
Kondo screening cloud, we have developed MC simula- 
tions. Since the nonlinear terms in Eq. are local, we 
integrate out all fields away from x = 0. Under a path- 
integral representation, we can rewrite C{x) as an aver- 



age over new fields qj{T) = v27r0j(O, t), where j — c,s 
and r is the Euclidean time, and over the impurity spin 
field S{t) — 2Sz{t) — ±1. The <Jx,y,z operators have 
to be treated dynamically as well, but from Eq. (^) it 
follows that the corresponding field is constrained to be 
(r) — hS{t) with /i ~ ±1. We find the formal result 

C(x) = --^ cos[2kFx]Wc{x)Ws(x) D(x) . (12) 
27ra 

The functions Wj{x) describe an algebraic decay ~ 
2;-9j/2 on scales a; <C xt, followed by a crossover to an 
exponential decay, 



Wjix) = I'^sinh 

TT 



2'KX 



D{x) = -( ^S{t ^ 0) cos 



The impurity average is now contained in 
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where the average is taken using the action 

j=c,s ui 

■ ^iJ/47rvF?Y.\^\\Siu;)\' . 
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Frequency sums run over Matsubara frequencies, and 
qj{uj) and S{uj) are the Matsubara components of the re- 
spective fields. Discretizing Euclidean time into N slices, 
Tj = jAt with At — P/N, the part Sj becomes 



e"""" = J™ T\{tiS,+i,S,+i\e^p[~ATHj{T,)]\fiS,,S,) , 

(16) 

where Hj{t) is the last part (~ J) of the Hamiltonian 
(^), with -\/27r0c/s(O) being replaced by qc/s{T). The 
matrix elements can be evaluated in closed form, with 
Gz parametrized by ^Si with Si — S{Ti) — ±1. Since 
exp[— 5*^] is negative for certain impurity spin paths, our 
simulation method has to deal with the conventional sign 
problem [ po[ . Fortunately, the sign problem is moderate 
except near T = 0. 
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FIG. 1. Monte Carlo data for D{x) at = 1/2 and 
I3u)c ~ 100. Statistical errors are of the order 5%. Notice the 
logarithmic scales. 
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FIG. 2. Numerical results for the plateau value Do as a 
function of J for Qc = 1/2 and two different temperatures. 
Statistical errors are of the order of the symbol sizes, and 
dotted and dashed lines are guides to the eye only. 

In Fig. I, MC data for D{x) are shown for several J at 
Qc = 1/2. For small J, the function D{x) exhibits a power 
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law a;'' for a; <C xt, where 6 coincides with the RKKY law. 
Far away from the impurity, D(x) reaches a plateau value 
Do in agreement with Eq. (^). Therefore the RKKY law 
is fully reproduced by our simulations. For large J, the 
numerical results display a different behavior. The func- 
tion D{x) now decreases to a small plateau value, and 
the RKKY law breaks down even at short lengthscales. 
From our numerical data, one has a complete breakdown 
of RKKY for J > J* with J* /2ttvf ~ 0.1. Furthermore, 
the numerical simulations predict the asymptotic expo- 
nent (1 -I- (7c)/2, since D{x) generally reaches its plateau 
value Do for x < xt- From Eqs. ( p^ ) and (p^, one then 
infers the asymptotic form of C{x), 

C{x) - cos[2kFx]{x/a)-'-^+s-^/^ , vf/Tk ^ x <t: xt , 

(17) 

which we have also verified by using lower simulation 
temperatures than in Fig. |l|. 

In view of Fig. |l|, it seems convenient to discuss the 
suppression of RKKY oscillations by Kondo screening in 
terms of Do. Numerical results for the plateau value Do 
at = 1/2 are shown in Fig. Taking some fixed 
J < J* and then going to low temperatures leads to 
an increase in Do. On the other hand, for J > J*, we 
observe a decrease in Do with lower temperatures. This 
indicates a crossover from a regime J < J* , where RKKY 
behavior is observed, to a non-RKKY regime J > J*. 
Finally, for the special value J — 2ttvf (Toulouse limit), 
i.e., J = 0, one finds the exact result C{x) = implying 
that Do ^ as J approaches the Toulouse limit. The 
correlation function C{x) vanishes identically since the 
Hamiltonian (|^) stays invariant under the transformation 
9s{x) —9s{x), whereas Eq. (||) changes sign. 

It is instructive to compare the asymptotic behavior 
( |l7| ) of C{x) with the Friedel oscillation of the charge 
density. Renormalization group and conformal field the- 
ory imply that in the strong-coupling limit, Sz and Sz(0) 
form a local singlet |l^ . This singlet decouples from the 
system and simply acts as an elastic potential scatterer 
in the unitary limit. The Friedel oscillation for that case 
is given in Ref. [Q. In a magnetic field B one obtains 
for spin a = ± 

P.{x) = ^ - , (18) 

where fc^ = kp +aB/AvF and — l/2gckp. Clearly, 
the Friedel oscillation and C{x) are both characterized 
by the same asymptotic exponent. 

As demonstrated in Ref. ||l^] , the Friedel oscillation can 
also be employed to determine the T = local suscepti- 
bility (|7|). This quantity is experimentally accessible in 
terms of the Knight shift. Using xi^) = d{sz{x)) / BE and 
{sz(x)) = X^cr ''^'°<t(^)/'^' obtain the leading asymp- 
totic behavior. 



A'KVFxi.x) = -(x/ao)^^"^'^)/^ cos[2kFx\ + 1 (19) 
_ 9c{l~9c) sin(2fc^^)(^/ap)-(i+gc)/2 ^ 

where ao — l/2gckF- Remarkably, for correlated conduc- 
tion electrons, the Knight shift actually increases with 
distance. A related behavior has been reported for a 
non-magnetic impurity in a Heisenberg chain pl| . 

To conclude, for correlated electrons, the RKKY inter- 
action exhibits only a slow algebraic decay. This implies 
that the usual logarithmic 2fci?-singularity of the ID sus- 
ceptibility is turned into an algebraic divergence. Fur- 
thermore, there is an interesting crossover from RKKY 
to Kondo screening cloud behavior. Both are character- 
ized by different exponents, and both lead to a slower 
decay than in the noninteracting case. 
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